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Abstract: We consider the large N limit of four dimensional SU(N) Yang- Mills field 
coupled to adjoint fermions on a single site lattice. We use perturbative techniques to show 
that the Zfj center-symmetries are broken with naive fermions but they are not broken 
with overlap fermions. We use numerical techniques to support this result. Furthermore, 
we present evidence for a non-zero chiral condensate for one and two Majorana flavors at 
one value of the lattice gauge coupling. 
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1. Introduction 

Continuum reduction [1] holds in the 't Hooft limit of large N gauge theories in d > 2 
enabling one to obtain physical results in the infinite volume limit and zero temperature 
by working on a l d torus with I of the order of one or two fermi. Reducing I below a certain 
critical value will force the theory to go into the deconfined phase. This in contrast to the 
theory in d = 2 where one can set I = and work on a single site lattice. This is referred 
to as Eguchi-Kawai reduction [2] and it works in d = 2 since the physical theory is in the 
confined phase for all temperatures. 

A weak coupling analysis [3] on a single site lattice shows that the large N limit of 
SU(N) Yang-Mills theory breaks all Zfj symmetries associated with the Polyakov loops 
if d > 2. Fermions in the fundamental representation will not affect this argument but 
fermions in the adjoint representation can modify the results of the weak coupling analysis. 
A continuum analysis of the theory with adjoint fermions on R 3 xS 1 with periodic boundary 
conditions for fermions in the compact direction shows that the Zjq symmetry is not broken 
in that direction [4] . An analysis on S" 3 x S 1 also shows a region where the Zn symmetry 
is not broken [5] . A lattice analysis of the same theory with Wilson fermions indicates 



- 1 - 



that one can reduce the compact direction to a single site on the lattice and still maintain 
the Zjv symmetry [6, 7, 8]. 1 

The natural question that follows is if the Zfj symmetries remain unbroken on a single 
site lattice in d > 2 dimensions. Recent numerical analysis with Wilson fermions in d = 4 
has shown this is most likely the case for a wide range of quark masses [10]. 

In this paper we address the same question as above in d = 4 using naive fermions and 
overlap fermions. The reason to study naive fermions is the apparent absence of doublers 
on a single site lattice since we only have the zero momentum mode with periodic boundary 
conditions. But, we will show using weak coupling analysis and numerical analysis that 
naive fermions break the symmetries but overlap fermions do not. The reason for this 
difference is a certain Z<i symmetry that is broken by naive fermions. 

Fermions play a dynamical role and we opted to use the hybrid Monte Carlo algorithm 
for generating gauge field configurations. Since we are on a single site lattice, the size of 
the fermion matrix is 4(iV 2 — 1) x 4(iV 2 — 1) and we compute the fermionic force exactly. 
We perform a full inversion of the overlap Dirac matrix using an exact diagonalization of 
its kernel, the Wilson-Dirac matrix. We present the details of the underlying hybrid Monte 
Carlo algorithm whose computational cost scales like iV 6 . We will show in this paper 
that one can use this algorithm to compute physics quantities like the chiral condensate. 
Detailed physics results will be presented in [11]. 

2. The single site model 

The partition function on a single site lattice is given by 

Z = J [dU]e~ s = J [dU]e- s °- s f, (2.1) 

where 

4 

S g = -bN Tr [u„U u Ulut\ , (2.2) 

is the standard Wilson gauge action. With / flavors of Dirac fermions, the contribution to 
the action from fermions is 

S n ,o = -/ log det (2.3) 

where "y^H nfi is the naive (n) or overlap (o) Dirac operator. The N x N matrices U^; 
/i = 1,2,3,4 belong to SU(N). The lattice gauge coupling constant is b = The 
fermions are in the adjoint representation and they couple to the gauge fields, V^, given by 



yab = ZTl 

n 2 



1. 



T a U„T h Ut 



(2.4) 



where T a , a = 1, • • • , (iV 2 — 1), are traceless hermitian matrices that generate the su(N) 
lie algebra. We have chosen the generators, T a , to satisfy 



TlT a T b = 2§ ab^ ( 25 ) 



1 There is a related paper on this subject [9] and the reader is referred to [7, 8] for a clarification of 
possible discrepancies. 
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The fully anti-symmetric structure constants are defined using 



[T a ,T b ] = J2iffT c . (2.6) 

c 

Both H n and H Q are 4(iV 2 — 1) x A(N 2 — 1) hermitian matrices and correspond to 
naive Dirac fermions and overlap Dirac fermions respectively. The determinant of H n o is 
positive definite and therefore the logarithm is well defined. Furthermore, we will show in 
the following two subsections that there exists a hermitian matrix S such that 

EiJ nj0 E = H* a , (2.7) 

which implies that all eigenvalues of H n ^ a are doubly degenerate reflecting the adjoint 
nature of the fermions. In addition, both naive and overlap fermions obey chiral symmetry 
and therefore the eigenvalues of H nfi will come in ± pairs. Therefore the factor, /, in front 
of can be an integer (single Dirac flavor) or half-integer (single Majorana flavor) for all 
values of N ? Since / appears only as a multiplicative factor in the fermion action, we may 
take the number of fermions flavors, /, to be any real number and we will do so in this 
paper. 



2.1 Naive fermions 

The matrix H n is given by 



where [i is the fermion mass. The chiral matrix, C, is 

C=lY,^(V»-Vt); (2.9) 

-a;)— a G-j— a ■)• <-> 

Note that 

o"2^<72 = -<V (2-11) 
Since are real matrices, it follows from (2.11) that 

C* = -a 2 Ca 2 . (2.12) 

If we define 

then (2.7) follows from (2.12). 

2 Note that one should have written 7s-ff n ,o in (2.3) but this is the same as writing H n . as long as / is 
an integer multiple of |. 
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2.2 Overlap fermions 

The hermitian Wilson Dirac operator is given by 

= (4 - m) 75 - £ + «,t Vj) (2.14) 

with m being the Wilson mass parameter and 

1/1 —o"„ \ 

(2.15) 



2v4 -i 



Using the definition of £ from (2.13), we have 

Su^S = w*; S75S = 75, (2.16) 
and therefore it follows from (2.14) that 

SifS = H*. (2.17) 
The hermitian massive overlap Dirac matrix, Hq, is defined by 

# = ~[(l + //)75 + (l-/i)e(ff)], (2.18) 
where \x £ [0, 1] is the bare mass. (2.7) for overlap fermions follows from (2.17). 
2.3 Fermion boundary conditions 

We have assumed periodic boundary conditions for fermions in the previous two subsec- 
tions. Other choices of boundary conditions that do not generate a U(l) anomaly amount 
to replacing by V^e 1 n with integer valued [12]. Physical results are expected to 
depend on the choice of boundary conditions. This is in contrast to the case of large N 
gauge theories coupled to fundamental fermions. In that case, the e l n factor can be ab- 
sorbed by a change of gauge fields that only changes the Polyakov loop and not the action. 
If the Zjy symmetries are not broken as is the case in the confined phase, this change will 
not affect physical results. 

3. Weak coupling analysis 

Our aim is to study whether the Zfj symmetries are broken in the weak coupling limit. 
We follow [3] and perform the weak coupling analysis by decomposing according to 

£/ M = e ia "ZV~ mM ; D^ = e<Sij. (3.1) 

Keeping 6^ fixed, we expand in powers of a^. The lowest contribution to S g comes from 
the quadratic term in [3] and the lowest contribution to Sf comes from setting a At = 0. 
Each Vn has N ( N ~ l "> two by two blocks of the form 



(3.2) 



cos(^ - 91) sin(f£ - 61) 
-sin(0i-0£) coB(0*-0fc 



- 4 - 



with 1 < i < j < N. The remaining (N — 1) x (TV — 1) matrix is a unit matrix. Therefore, 
the gauge field effectively has (N — 1) zero momentum modes and N(N — 1) non-zero 
momentum modes of the form e l ^~ e ^ with 1 < i ^ j < N. 

The computation of the fermion determinant reduces to a free field calculation at this 
order and the result is 



S n ,o = -4/ £ In A n , o (0* - & + <j>) - 4(iV - 1)/ In A n , o (0) 



(3.3) 



where e 1 ^, </>^ = is the phase associated with the boundary condition in the ^ 

direction. The eigenvalues, ±A(p), are two fold degenerate and the explicit expressions are 

\ n (p)= //x 2 + ^sin 2 ^ (3.4) 



for naive fermions and 



K(p) = 



l-/x 2 



2£,sin 2 f- 



2E M sin 2 f-^ + I>in 2 ^ 



(3.5) 



for overlap fermions. The complete result from fermions and gauge fields is 

1 



S = £ In 



4/lnA n , o (0 l -^ + , 



4(AT-l)/lnA n , (</>). (3.6) 



Independent of the actual values of 9* the fermion eigenvalues will have (N — 1) zero 
modes with periodic boundary conditons when the mass is set to zero. If all the 9^ are 
different for each fx, then the fermions should not have exact zero modes when we set 
equal to (o 1 ^ — 9j^j with i ^ j. If the fermion spectrum has more than (N — 1) zero modes, 
we will refer to the extra ones as doubler zero modes. 

In order to find the minimum of S, we consider the Hamiltonian 



(3.7) 



For large f3, the Boltzmann measure e~ H will be dominated by the minimum of S. We can 
perform a HMC update of the ir, 9 system to find this minimum. The equations of motion 
for our H are 

(3.8) 
(3.9) 



and 



dS 
d9j, 



-p {fi - 2/f; d ) . 



The gauge contribution to the force is 

*S = E 



sin 



4) 



sin 



2 h 9i-9i 



(3.10) 
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The naive fermion contribution to the force is 



F' 

n 



EE 



sm 



a=± j^i V + °a 

and the overlap fermion contribution to the force is 



EE 



11 11 1 • 
SmpjiaCa ~ 2 sm 



2?; 



i+M 2 



where 



^sin 2 ^ ± 



It is clear from (3.9-3.12) that 



2 ?v± 



sm 



+ c l l 



m 



If we start with 



i 

then it will remain zero. Then it follows from (3.8) that 



If we start with 



£«4 = °. 



then it will remain zero. These are just the conditions for remaining in SU(N). 
A choice for the order parameters associated with the Zfj symmetries is [3] 



AT2 



(3.11) 



(3.12) 



(3.13) 
(3.14) 

(3.15) 

(3.16) 
(3.17) 



(3.18) 



If = 2 5 then the symmetry in that direction is not broken. If are uniformly 
distributed in a width a < 2tt, then 



lim P u 



1 



2 a 

— sm — 
a 2 



(3.19) 
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3.1 Naive fermions break the Zfj symmetries 

We assume periodic boundary conditions and set <f>n = in (3.6). We pick one value of 
iV and (3 and calculate P^ as a function of /. In order to clearly see symmetry breaking 
we use rotational symmetry on the lattice and choose to label our directions such that 
P\ < P2 < P3 < -P4 for each configuration in our thermalized ensemble. We set /U = 0.01 
to avoid potential singularities that could occur for the massless case. The plot for three 
different choices of ((3,N) is shown in Fig. 1. The value is well below \ and it seems to 
approach a value for large / that is consistent with a = tt in (3.19). Since all three choices 
for (/3, N) give consistent values, we have plotted the results for (3 = A and N = 23 in 
Fig. 2. The small deviations one sees between the four different P^ are simply due to our 
ordering scheme combined with finite N effects. 

With N = 23, we expect 22 exact zero modes for \{p) as explained in section 3. We 
plot the average of the smallest eigenvalue, \{p) with p^ = (j) 1 ^ — 9j^j and i ^ j, in Fig. 3. 
It is clear that this eigenvalue is non-zero for all values of / indicating that there are no 
doubler zero modes. The reason for the breaking of the Zfj symmetries can be understood 
by looking at the total action obtained from (3.4) and (3.6): 



¥=3 



sin 



2/J> 

¥3 



sin 2 (^-^) 



(3.20) 



The fermionic contribution cannot separate 6^ = 6^ from 6^ = 9^ + tt implying that the 
fermion contribution alone will result in a distribution of eigenvalues restricted to a width 
of TT. 



3.2 Overlap fermions do not break the symmetries 

Contrary to naive fermions, overlap fermions do separate 6^ = Oft from 6^ = 0^ + tt as 
is evident from the presence of the Wilson factors in (3.5). Therefore, we do not expect 
overlap fermions to break the Zfj symmetries. A plot of P^ for several values of / at 
N = 23, (3 = 1, and m = 5 with /i = 0.01 in Fig. 4 shows this to be the case. The very 
small deviation close to / = \ is a consequence of finite N effects. Since we do not expect 
to go beyond iV = 23 in the full simulation of the model, this plot will serve as a guide to 
what one can expect in a full simulation. 

Since there are no doubler zero modes, we do not have any restriction on the values 
for the Wilson mass, m, used in the Wilson-Dirac kernel as described in section 2.2. But, 
we cannot make it arbitrarily large since one can see by a direct computation that the 
large m limit of overlap fermions is naive fermions [13]. A plot of P^ as a function of m is 
shown for j = \ and / = 1 in Fig. 5. It indicates that 3 < m < 8 is an appropriate range 
of values of m where the Zfj symmetries are not broken for / = 3 and that range only 
gets bigger as / increases. A plot of the lowest positive eigenvalue of H w in Fig. 6 and the 
smallest eigenvalue, X Q (p) with p^ = ^6^ — Oj^j, i 7^ j, in Fig. 7 shows that there are no 
doubler zero modes in this range of m. We can use this range of m for our full numerical 
simulation with overlap fermions. One cannot take the Wilson mass parameter to zero in 
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Figure 1: Naive fermions: Plot of as a function of / for three different choices of (f3,N) and 
/i = 0.01. 

the weak coupling limit. This is intimately tied to the fact that the Zjy symmetries are 
not broken and ^ 1 in the weak coupling limit. 

3.3 Effect of boundary conditions 

In order to study the effect of boundary conditions on symmetry breaking, we focus 
on / = g and / = 1. We set 4>^ = for \i = 2,3,4 and varied (\>\ by setting it to equal 
to with k an integer in the range < k < N/2. The plot of as a function of tpi is 
shown in Fig. 8. We see that the symmetry in the /i = 1 direction is broken if (pi > ^. 
This seems to be the case in the limit of large N and seems to be roughly independent of 
/. Furthermore, the Zn symmetries in the other three directions with periodic boundary 
conditions are not broken. This result could help us force feed momentum for quarks in 
the adjoint representation. In order to pursue this, we need to study the effect of (pi on the 
chiral condensate and see if chiral symmetry is restored when the Zn symmetry is broken 
and if the chiral condensate is independent of the value of (pi when the Zn symmetry is 
not broken. 

4. Hybrid Monte Carlo algorithm 

One could use the standard Metropolis algorithm to study the single site model but the 
computational cost would scale like iV 8 [10]. We will use the Hybrid Monte Carlo algo- 
rithm [14] to numerically study the model on a single site lattice. Fermion matrix inversions 
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Figure 2: Naive fermions: Plot of P M as a function of / at iV = 23, /3 = 4, and fi — 0.01. 



dominate the computational cost of this algorithm and will therefore scale like iV 6 . We 
set Hp as momentum variables conjugate to lf$ with the condition that H^ is a traceless 
hermitian matrix. Therefore, we can write H^ as 



(4.1) 



The Hamiltonian for the Hybrid Monte Carlo algorithm is 



(4.2) 



Since the algorithm has not been used in the past for gauge theories with adjoint fermions 
on a single site lattice and since we do not use the conventional pseudo-fermion algorithm 
for dealing with the fermionic determinant, we present the necessary details in this section. 
The equations of motion for £/„ are 



Setting = results in 



dUp 
dr 



dH, 



iH^. 



dSn 



(4.3) 



j, i dSn.i 
_| » _| - 

dr dr 



(4.4) 
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Figure 3: Nai've fermions: Plot of smallest eigenvalue, X(p) with = (6' 1 — 8fy and i ^ j, for 
N = 23, = 4, and fi = 0.01. 



where we can write 



' As '"" ^Tr //„/;;;•" . (4.5) 



A simple calculation results in 

Dl = -zbN E [t^t^C/jL^ + UJJ$UlU v - UtU^Ul - UMMK} ( 4 - 6 ) 

V 

The derivative of is 

-^- = -£ Hpcfrv* = ~J2 s; d v^ (4.7) 

cd d 

where 

H a ^ d = J2 H ^fd C , (4-8) 

c 

and we have used (2.6) and (4.1) . Using the anti-symmetry of the structure constants, it 
follows that are anti-symmetric real matrices. In analogy with (4.3), we can write 



- -HM. (4.9) 
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Figure 4: Overlap fermions: Plot of as a function of / at N = 23, (3 — 1, m = 5, and = 0.01. 



4.1 Derivation of D™ for naive fermions: 



We start by noting that 



dC , dC}_\ 
m 2 + C& \dr dr 



and 



Let us define 



dC _ _1 



Mp, = -tr 



m 2 + CC^ av 



; m\ v = m vii , 



(4.10) 
(4.11) 
(4.12) 



where tr stands for trace only over the spin index. Using (2.11) and (2.12) we see that 
My V are real matrices. We can rewrite (4.10) as 



dSn 

dr 



where 



Ay = [V^B^ + BfjV*] - transpose, 
are real anti-symmetric matrices and 

b, = ]T(K - v$)m m 



L fJLV 



(4.13) 
(4.14) 
(4.15) 
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Figure 5: Overlap fermions: Plot of as a function of m at N = 23, (3 = 1, and /i = 0.01 for two 
different values of /. 

This brings us down to the form we need, namely, 

f^ = £W^ (4.16) 
n 

where 

Dl = i f -Y,Al b [T*,T h ). (4.17) 

ab 

4.2 Derivation of D° for overlap fermions: 

Let the eigenvalues of the Wilson-Dirac operator, H, given in (2.14) be 

H\Xi) = \i\Xi); \i>0 
H\ui) = 0Ji\uJi); tOi < 
{Xi\Xj) = Sij; (uJi\uj) = 5if (Xi\ujj) = 0. (4.18) 

The dimension of H is 4(iV 2 — 1) x 4(iV 2 — 1) and we will assume that there are an equal 
number of positive and negative eigenvalues. 3 If follows from (2.17) that all eigenvalues 
of H are doubly degenerate since the pair of eigenvectors (E|Aj)*, \Xi)) have the same 
eigenvalue and the same is the case for the pair (T,\coi)* , \cui)). 
The fermion action is defined by 

S = -fTilogH , (4.19) 
3 This is a reasonable assumption since we can restrict ourselves to a fixed topological sector. 
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Figure 6: Overlap fermions: Plot of the lowest positive eigenvalue of H as a function of m at 
N = 23, /3 = 1, and /x = 0.01 for two different values of /. 



and 



It follows from (2.18) that 



We can write 



dSo = 1 dH 

dr H dr 



dH _ 1 - jx de{H) 
~dV ~ ~~2 dr~ ' 

e (F) = ^|A i )(A i |-^K)(a; i |. 



Therefore, 

de(H) 
dr 

and 



\- d\\j) ^ d(\j\ \- d\uj) \ - 



d|A 4 ) _^(A i |^|\) 



The above equations imply 



E 



E 



A, -A, 
u>i — Xj 



|A;> + £ 



|A;> + £ 



Aj — 

a;, — a;,- 



■J/- 



de{H) 
dr 



|A fc ) =2^ 



Afc — Uj 



(4.20) 

(4.21) 
(4.22) 

(4.23) 

(4.24) 
(4.25) 
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Figure 7: Overlap fermions: Plot of smallest eigenvalue, X(p) with — (O 1 ^ — 6fy and i ^ j, for 
N = 23, /3 = 1, and /i = 0.01. 



de(#) , 
dr 



— A,- - uj k 



(4.26) 



Inserting (4.21) into (4.20) and using the above equations, we get 



dS 



/(l-M)Tr 



dH 



G + G f ) 



with 



3,k 



\*j)(*j\-H;\ u k)(wk\ 



idfr — Xj 

The symmetries in (2.17) and (2.7) imply that 

SGS = G*. 

Equations (4.27), (2.14) and (4.9) can be combined to give 
dS 



dr 



/(l - rfTrHn 



vJG + Gtjup-uUG + Cftjv* 



Let 



= tr (G + Gt)w M , 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



(4.31) 
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Figure 8: Overlap fermions: Plot of P n as a function of 0i for three different values of (/, N) with 
@= 1 and n = 0.01. 

where tr stands for sum over the spin indices. Using (2.16) and (4.29) we can see that 
are real. With the definition of B^, (4.30) becomes 

A Q 

= /(l - riTrHp [V^ - B l X] = /(I " ^TrS^, (4.32) 

where 

-V = VfiBfj, - B l X (4-33) 

are real anti-symmetric matrices. The rest of the steps are identical to that of naive 
fermions. 

5. Numerical results 

In this section we present numerical results from our Hybrid Monte Carlo algorithm and 
show that they support the conclusions reached in section 3. 

Fig. 9 shows a plot of P^ as a function of / obtained using naive fermions. We have set 
b = 5, N = 11 and fi = 0.01. Naive fermions break the Zfj symmetries as expected from 
the weak coupling analysis. The values of Pa in Fig. 9 are consistent with the distribution 
of 9^ uniform with a width less than ir and this is evident for / = 1 in Fig. 10. 

Turning our attention to overlap fermions, we plot P^ as a function of m in Fig. 11 for 
6 = 7, N = 11, / = ^ and ji = 0.01. As expected the Zfj symmetries are not broken for a 
wide range of m. A plot of the lowest eigenvalue of H as a function of m in Fig. 12 shows 
that there are no doubler modes in the same wide range of m. 
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Figure 9: Naive fermions: Plot of P M as a function of / for TV = 11, b = 5 and = 0.01. 



If chiral symmetry is broken, we expect the low lying eigenvalues of the overlap Dirac 
matrix to scale like jmzzr- Furthermore, we expect agreement with chiral random matrix 
theory predictions for the distribution of the low lying eigenvalues in accordance with the 
symplectic ensemble [15]. Instead of using the analytical formulas in [15] we found in 
convenient to numerically generate the eigenvalue distribution governed by the symplectic 
ensemble using Hybrid Monte Carlo techniques and compare with the numerical results 
obtained for this model. The chiral Random Matrix theory ensemble for a symplectic 
matrix, C = Yju 



is 



Z 



[dC^e 



' 5Z M [C/i 



[det H 



rmtj 



l/. 



H 



rmt 



/' 

at 



c 



(5.1) 



(5.2) 



with being a real square matrix. We expect one scale factor to relate the distributions 
of the eigenvalues of H Tmt and H Q . Since we want to eliminate the scale set by the chiral 
condensate, we focus on 

M, 

where < Ai < A2 are the first two non-degenerate eigenvalues. Fig. 13 shows a plot of r 
as a function of / obtained at = 11 and 6 = 5 with m = 5 and [i = 0.01. The results 
are also compared with the numbers one gets from chiral random matrix theory. There is 
agreement with chiral random matrix theory for f = \ and / = 1 but not for / = 2 or 
/ = 3. The case of / = | corresponds to the large N limit of super Yang-Mills and this is 
one case where chiral symmetry assures supersymmetry. 
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Figure 10: Naive fermions: Distribution of 9 l for N = 11, 6 = 5 and /j, = 0.01 at / = 2. 

One can try to understand the above behavior by considering the running of the 
coupling constant, a = j-jt with the scale set by the lattice spacing, a [18]. Using the 
two-loop beta function [16, 17], we have 

da /ll 4 \ , /34 32 A o . . 

where the first two coefficients listed are independent of the the renormalization scheme. 
Asymptotic freedom is lost when / > -j- . The two-loop expression predicts an infrared fixed 
point [19] if ^| < / < if with the fixed point coupling becoming smaller as / increases in 
this range. Our current data (Fig. 13) is consistent with the presence of a chiral condensate 
in the absence of an infra-red fixed point. One should keep in mind that one has to take 
the large N limit at a fixed b and then take the continuum limit, b — > oo, in order to arrive 
at a physically relevant conclusion. 

6. Future outlook 

We have provided compelling arguments for the preservation of Zfj symmetries in the large 
iV limit of Yang-Mills theories coupled to adjoint overlap fermions on a single site lattice. 
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Figure 11: Overlap fermions: Plot of P M as a function of m for = 11, b = 7, / = | and /i = 0.01. 

We have provided all the necessary details to perform a Hybrid Monte Carlo algorithm and 
we have shown numerical feasibility. We plan to study several different values of N in the 
range of 11 to 18 and we also plan to study several different couplings in the range of 0.35 to 
5.0 [11]. One can use a physical observable like the Wilson loop to study the beta function 
on the lattice and investigate the presence/absence of an infrared fixed point. In this 
context, it would be useful to perform a weak coupling expansion of the beta function on 
the single site lattice and confirm that one obtains the same result as on an infinite lattice. 
Perturbative computations on a single site lattice do not have the standard divergences and 
a perturbative computation of the Wilson loop will clarify the role of perimeter divergence 
on a single site lattice. Finally, one should extend the numerical results to other choices of 
boundary conditions that are not periodic. 
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